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ABSTRACT 

We present a 0.5-2.2/xm scattered light spectrum of the circumstellar disk 
around TW 7 Hya from a combination of spatially resolved HST STIS spectroscopy 
and NICMOS coronagraphic images of the disk. We investigate the morphology 
at the disk at distances > 40 AU over this wide range of wavelengths. We measure 
the surface brightness, azimuthal symmetry, and spectral character of the disk as 
a function of radius. We find that the scattering efficiency of the dust is largely 
neutral to blue over the observed wavelengths. We find a good fit to the data 
over a wide range of distances from the star if we use a model disk with a partial 
gap of 30% depth at 80 AU and with steep disk truncation exterior to 100 AU. If 
the gap is caused by a planetary companion in the process of accreting disk gas, 
it must be less than 20 M®. 
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1. Introduction 


TW Hva is the closest example of a star with a young, gas-rich protoplanetary disk. As 
such, it is an ideal target in which to study how disk evolution is coupled to planet formation 
because the disk can be studied at high spatial resolution and with great sensitivity. First 
detected from its infrared excess (Rucinski & Krautter 1983), the disk emits strongly at 
wavelengths longer than 3 /rm with a peak in its spectral energy distribution between 10 and 
100 fim (Weinberger et al. 2002). 

TW Hya may not be typical for its age. It is part of the eponymously named TW 
Hydrae Association, a collection of about two dozen stars, amongst which it sports the most 
massive, gas-rich disk of all. TW Hya is 54 ± 7 pc away in the new van Leeuwen (2007) 
Hipparcos catalog. The age of the TW Hya association is found from a variety of studies of 
the ensemble of stars, including measurements of the dynamics (8.3 ± 0.3 Mvr) (Weinberger 
2011; de la Reza et al. 2006), lithium depletion boundary (12 ± 8 Myr) (Mentuch et al. 
2008), and these in combination with pre-main sequence tracks (lOl^ 0 Myr) (Navascues 
2006). Associations just a bit older than TW Hya, such as /3 Pic, have no optically thick, 
massive protoplanetary disks. 

Despite reservations about its representativeness, intensive studies at every wavelength 
and spectral resolution have been used to try to understand the TW Hya disk’s structure 
and composition. We summarize some of the main findings here. Resolved emission at 7mm 
indicated the central disk is mostly clear inside of 4 AU (Hughes et al. 2007). However, 2 /mi 
emission resolved by the Keck Interferometer show that an optically thin dust disk comprised 
of small grains comes to within 0.06 AU of the star. The spectral energy distribution, 
particularly the large amount of emission at submm through cm wavelengths (Weintraub 
et al. 1989; Wilner et al. 2000: Weinberger et al. 2002; Wilner et al. 2005; Andrews et al. 
2011) indicates substantial grain growth to at least cm sizes. Resolved CO line maps show 
the disk in Keplerian rotation with an inclination of 7 ± 1° beyond 50 AU (Qi et al. 2004; 
Andrews et al. 2011). 

The disk has been spatially resolved in scattered light in the visible and near-infrared 
(Krist et al. 2000: Weinberger et al. 2002; Apai et al. 2004; Roberge et al. 2005). These show 
that the optically thick disk extends to at least 280 AU. They also show asymmetries and 
changes in surface brightness of the disk in the inner 150 AU. The presence of both small 
and large grains throughout the 4 AU - 200 AU disk and the presence of small grains in 
a region where they should quickly be removed, suggests that they are being regenerated 
through collisions at multiple locations. 

The disk chemistry has also been probed in a spatially resolved manner in submm lines 
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(Qi et al. 2006). For example, there appears to be ongoing deuterium fractionation in the 
outer disk, and that suggests pristine nebular material is not preserved, as is often assumed 
for comets. 

We have undertaken a multiwavelength visible to near-infrared study of scattered light 
from the disk in order to address several questions of interest to planet formation models. 
Observations in the visible to near-infrared can detect ices and silicates due to their broad 
absorption features and organics due to their red slopes. These types of observations are 
routinely applied to comets and Kuiper Belt objects, which are thought to be the planetesimal 
remnants of ci reruns tellar disks. In addition, the spectral scattering efficiency of the disk 
can constrain the grain sizes at the disk surface. The mixture of grains should reflect a 
combination of vertical mixing from the midplane, where large grains are presumably formed, 
radial transport, and collisions. Finally but by no means least, we wish to understand the 
vertical structure of the disk and whether it shows evidence for forming planets. 


2. Observations 

We took coronagraphic images with the F171M, F180M, F204M, and F222M filters 
(central A=1.71, 1.80, 2.04, and 2.22/nn respectively) on 09 May 2005 with the NICMOS 
camera 2 for TW Hya and the PSF reference star, CD-43°2742 as part of Program GO 10167. 
The observations include direct images of both stars outside the coronagraphic hole with 
short exposures for point source photometry, as well as longer exposures for coronagraphic 
high contrast imaging. The instrumentally calibrated and reduced images discussed in this 
paper were created from the raw NICMOS multiaccum exposures following the processing 
methodolgy described by §3 of Schneider et al. (2005) and references therein. 

For photometric analysis, each calibrated direct image was used to determine the total 
photometry of the star and empirically determine the scaling ratios between TW Hya and the 
PSF reference in each filter band. The three images for each star in each filter were located 
at different positions on the detector. We used a median combination of the three dither 
points to create a final image of each star to derive a ratio for scaling and for photometry of 
TW Hya. We used a 16.5 pixel radius circular aperture to determine the photometry. The 
background in the images is zero, so no background annulus was used. The individual dither 
points were used to get a rough estimate of the uncertainty in the ratios and photometry. 
Table 1 lists the photometry of TW Hya in each band with uncertainties and the scaling 
ratios for each filter. 


In order t.o determine the best subtraction we minimized a chi-squared metric on a 



region of the target image dominated by the star’s diffraction spikes. We assumed that 
good subtraction of the diffraction spikes corresponded to the best subtraction of the PSF 
within the region of interest (Schneider et al. 2001 ) . We iteratively created subtractions for 
combinations of scaling and pixel offsets until we found an image that produced the lowest 
chi-squared measure. We searched within 1-cr of the scaling ratios we determined above and 
within ±1 pixel to find the best x and y pixel offsets. 

To quantify the systematic effects on the photometry, we repeated the subtractions 
varying the PSF scalings and offsets by ±f a from the minimum chi-square solution found 
above. Using a circular photometric aperture matched to the size of the disk, we found 
the standard deviation in the disk flux densities from this suite of subtractions. We then 
propagated this uncertainty into the total uncertainty in the flux density of the disk per 
pixel. 

We observed TW Hya in each medium band filter at two distinct spacecraft orientations 
or roll angles. This is essentially an azimuthal dither that allows true objects within the field 
of view to be distinguished with instrumental artifacts that do not rotate with a change in 
orientation. The reference star, CD-43°2742, was subtracted from both roll-angle images to 
create two separate subtraction images for each filter. The images were then geometrically 
corrected for the slight optical distortion of the NICMOS camera 2 at the coronagraphic 
focus. We used the x-direction pixel scale of (T.'07595/pixel and the y-direction pixel scale 
of Cf/07542/pixel to create an image with pixels that have the y-direction plate scale in 
both directions. The geometrically corrected images were rotated to a common celestial 
orientation using the rotation centers given by the flight software in the raw data file headers. 
Figure 1 shows the resulting PSF subtracted images of the TW Hya disk taken in the medium 
band filters. 

Additional observations of TW Hya with NICMOS were performed as a part of programs 
7226 and 7223 with HST (Weinberger et al. 2002). The observations in F110W and F160W 
were recovered from the archive, reduced in the same manner as the medium band data, 
and median combined. Archival PSF reference stars for the images were subtracted from 
the target observations. The F110W' reference was r 1 Eridani, while the F160W reference 
was G1 879. We followed the procedure of Weinberger et al. (2002) for scaling the PSFs, but 
followed the above procedure for subtraction. 

TW Hya has been observed with STIS using a combination of wedge coronagraphic 
images as well as spatially resolved spectroscopy in the visible. The reduction for the data 
is detailed in Roberge et al. (2005). We used the same datasets that were presented in that 
work, with one exception. The spatially resolved STIS spectroscopy was re-extracted in O'.' 35 
distance bins to match our photometric apertures. Figures 1 and 2 shows the STIS and 
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XI CM OS images of TW Hya. 


3. Results 

The observation of the disk in spatially resolved scattered light over a range of wave- 
lengths allows a powerful look into both the morphology of the disk and the scattering 
efficiency of the disk’s dust grains as a function of wavelength and distance from the star. 
In this Section we measure the surface brightness profiles, azimuthal asymmetries, and the 
scattering efficiency of the dust in the TW Hya disk. 


3.1. Radial Surface Brightness Profiles in the Medium Band Filters 

The radial structure of TW Hya belies its apparent smooth decrease in surface brightness 
as a function of radial distance. In fact TW Hya has some radial structure that deviates 
from what one would expect from a simple flared disk. We investigate the radial surface 
brightness profiles in our medium band NICMOS data and compare the behavior in these 
wavelengths to that seen in the shorter wavelength data. 

Figure 3 shows the F171M through F222M surface brightness profiles. In all of them 
we see the characteristic shift in behavior between ~80-130 AU that is seen in the visible, 

1.1, and 1.6/nn, namely a change in the slope of the surface brightness. This feature is now 
seen in all the wavelengths of light in which TW Hya is observed, indicating that this is a 
feature caused by morphology rather than composition. We interpret this feature as a gap 
or depression in the disk, which we explore further in §4.2. 


3.2. Azimuthal Surface Brightness Asymmetries 

Roberge et al. (2005) measured the surface brightness as a function of azimuth for the 
50CCD image of TW Ilya’s disk and found it possessed a significant asymmetry between 
65-140 AU consistent with the asymmetry caused by a disk inclined to the line of sight with 
forward scattering dust grains. The magnitude and position of the asymmetry corresponded 
to a measured PA of the brightness maximum of 233.6°±4° and an inferred inclination of 
15l|f under the assumption of a Henyey-Greenstein scattering phase function for the grains 
with an asymmetry parameter g of 0.5. They found no significant asymmetry for larger radii. 

In this paper we re-analyze the 50CCD images as well as the images at other passbands 



to search for similar asymmetries For each passband, we took the final subtracted images 
and constructed photometric apetures in the shape of concentric annuli as a function of 
increasing radius with equal spacings of 07357 (20 AU). Each annulus was then further equally 
subdivided into azimuthal bins with width 20°. For each image a total of 6 annuli were used, 
from 0"88 (50 AU) to 2 / /31 (130 AU). For the 50CCD image, there was enough signal-to- 
noise to extend out to 3791 (220.7 AU). We therefore had a measure of the surface brightness 
distribution between 50-130 AU for all passbands, and a measure of the distribution for the 
50CCD image and F110W image between 130 and 230 AU. 

Certain areas of each image were masked out due to the presence of residual diffraction 
spikes or in the case of the 50CCD image, the presence of the occulting wedge. For apertures 
where masked pixels comprised more than 2/3 of the aperture, no surface brightness was 
calculated. Each annulus was then scaled to the median of the brightest annulus to provide 
equal weight to the surface brightness in any one azimuthal bin at each radius. The final 
azimuthal surface brightness distribution was then constructed by taking a median of the 
first six annuli for all images and another distribution was constructed for the 50CCD and 
FI 10W image of the disk from the final four annuli. For the uncertainty at each point 
we took the calculated photometric uncertainty or the standard deviation of the azimuthal 
points, whichever was larger. 

Figure 4 shows the 50CCD and F110W results. Between 50 and 130 AU the asymmetry 
is similar in the two passbands. Beyond 130 AU, the character of the asymmetry is different: 
the minimum to maximum ratio of the asymmetry is smaller and the position of the maximum 
is not clearly aligned with that of the inner disk. The distance of 130 AU corresponds to a 
break in the radial suface brightness of the disk. For a comparison, we also overplot our best 
fitting scattered light, models of the disk, described in more detail in §4.4, for inclinations of 
7° and 15°. 

Figure 5 shows the F160W, F171M, F180M, and F222M azimuthal brightness profiles. 
We neglect the F204M data because it has lower signal-to-noise. For each of these four 
images, no significant azimuthal asymmetry is revealed. Whatever is the cause of the asym- 
metry at shorter wavelengths, it is not detected with significance in the longer wavelength 
observations. We discuss the potential causes of this in §6. 


3.3. Scattered Light Spectrum 


The scattered light spectrum of the disk around TW Ilya is a combination of the input 
stellar spectrum and the intrinsic spectrum of the disk. We removed the stellar spectrum of 
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TW Hva by dividing the measured surface brightness by the flux density of the star in each 
filter. We divided the STIS spectrum of the disk by the point source spectrum of the central 
star. 


TW Hya is too bright to be directly imaged in the F110W and F160W filters as well 
as the unfiltered 50CCD observation. We therefore estimated the direct photometry of TW 
Hya in the unfiltered 50CCD observation, as well as the F110W and F160W filters using the 
CALCPHOT task of the SYNPHOT 1 package in conjunction with Kurucz models tabulated 
in the SYNPHOT library. We used a solar metallieity, log g—4 model of a 3925 K star, 
assuming a luminosity of 0.25T0.06 L 0 based on a best fit match to TW Hya based on 
its estimated V magnitude during the STIS observation and parallax from the Hipparcos 
satellite and J, H, and K s photometry from 2MASS. We compared the predicted F171M, 
F180M, F204M, and F222M photometry with the SYNPHOT values-they agreed to within 
a, few percent of our measured photometry. 

Figure 6 shows the total reflectance spectrum of the disk from 0.5-2. 22pm and encom- 
passing radial distances of between 50-230 AU. This data includes the average over the same 
distances of the STIS spatially resolved spectra of TW Hya from Roberge et al. (2005), nor- 
malized to the STIS photometric point. The overall spectrum is relatively neutral between 
0.5-1. 6pm, becoming more blue at longer wavelengths. No significant signature of sharp ab- 
sorption lines are seen in any band. At longer wavelengths the disk becomes slightly bluer, 
possibly because of the presence of water ice in the observed dust. 


4. Modeling the TW Hya Scattered Light Disk 

In this Section we provide a model for the observed spectrum and morphology of the 
disk. This requires a model of the structure of the disk, discussed in §4.1. We infer that the 
depression we observe between 60-120 AU is caused by a “gap" in the disk. The origin of this 
feature will not be constrained by these observations, but a likely candidate is a planetary 
object that has succeeded in opening a gap. The limits to such a candidate are discussed in 
§5. However, other possibilities exist, which we outline in §6. 


1 http://www.stsci.edu/resources/softwareJiardware/stsdas 



4.1. Initial Disk Structure 


The disk models are generated as described in Jang-Condell (2008) and Jang-Condell 
k Turner (2011). The details of the radiative transfer modeling is also described in Jang- 
Condell k Sasselov (2003) and Jang-Condell k Sasselov (2004). The structure of the un- 
perturbed planet-less disk is generated in a two-step procedure: (1) we calculate a locally 
plane-parallel two-dimensional model for the entire disk, then (2) we take a radially limited 
slice of this disk and refine its structure, this time taking the full three-dimensional (3D) 
curvature of the disk into account. In both steps, we iteratively calculate the vertical density 
and temperature structure of the disk including radiative transfer and under the assump- 
tion of hydrostatic equilibrium. The main heating sources are stellar irradiation and viscous 
heating. 

In step (1), we use the same formalism developed by Calvet et al. (1991) and D'Alessio 
et al. (1998, 1999), with some simplifying assumptions. We calculate the disk from 0.25 to 
256 AU, with radial bins spaced by factors of \/2. At each radius, we assume that the disk 
is locally plane-parallel. 

In step (2), we remove the assumption of local plane parallelness and calculate radiative 
transfer in three dimensions as a numerical integration over the surface of the disk. Although 
we assume that the overall disk structure is axisymmetric, the calculation of the radiative 
transfer is done in 3D to include the curvature of the disk both in vertical height above the 
midplane (z) and in the azimuthal direction ( 0 ). 

We adopt for our stellar and disk parameters those used by Calvet et al. (2002): effective 
temperature of T* = 4000 K, radius i?* = 1 U . . and accretion rate of 5 x 10“ 10 M 0 yr _1 . This 
gives a total stellar luminosity of L — 0.23/. . . One exception is stellar mass. Instead of 0.6 
M 0 , as chosen by Calvet et al. (2002), we choose a mass of 0.4 M© The mass of the central 
star TW Hya is not well known: it is inferred from the spectral type and pre-main sequence 
evolutionary tracks. The canonical literature spectral type is K7 (Webb et al. 1999) but 
recent infrared spectroscopy concludes the star is M2. 5 and T e ff=3400 K (Vacca k Sandell 
2011). Most of the prior SED arid sub-mm spectral cube modelling has assumed M*=0.6 
M© (Calvet et al. 2002; Qi et al. 2008), but the new spectral type implies a stellar mass 
of 0.4M©. The main affect on our models of reducing the stellar mass is to make the disk 
puffier, hence able to intercept more stellar light, and therefore brighter, bringing it more 
into agreement with the measurements. We assume that the dust is well-mixed with the 
gas and constant throughout the disk, with a dust-to-gas ratio of 0.0138 and a grain-size 
distribution of n(a) oc a~ 3 5 . 

Calvet et al. (2002) found that large maximum grain sizes ( a max ) give good fits to the 



SED of TW Hya at long wavelengths (A > 1 cm). Since there is a degeneracy between 
maximum grain size and viscosity parameter ( a ) in their models, they found that good fits 
to the SED of TW Hya are obtained with a max = 1 mm, 1 cm, or 10 cm, with a = 5 x 10 — 4 , 
3 x 10~ 4 , or 1 x 10'~ 4 , respectively. This is because large grains are effectively invisible 
to optical and infrared wavelengths. The a parameter must be adjusted as the maximum 
grain size changes, because the overall mass of the disk increases with increasing grain 
size, assuming a fixed gas-to-dust ratio; however, the stellar mass accretion rate must be 
kept fixed. We find a similar degeneracy in our models, so we simplify our calculations by 
adopting a max = 1 cm and a — 2.5 x 10 ” 4 . 

On the other hand, the scattered-light appearance of the disk is sensitive to the minimum 
grain size assumed, a min , and the composition of the grains. The wavelength dependent 
opacities are calculated using a Mie scattering code, miex, that includes distributions of 
particle sizes and can account for large dust particles (Wolf & Voshchinnikov 2004). For a 
model with a given minimum grain size, we recalculate the entire disk model from the initial 
conditions in order to be completely self-consistent. 


4.2. Gap and Truncation 

The gap in the disk is modeled as an ad hoc axisymmetric density perturbation param- 
eterized by a Gaussian with adjustable width w and depth d centered at 80 AU. To model 
the truncation of the disk, we introduce an exponential cutoff at a knee k. as appears in the 
self-similar solution for an accretion disk as derived in Hartmann et al. (1998). If E 0 (r) is 
the unperturbed disk surface density, then the new surface density profile is 

E(r) = E 0 (?’){1 — dexp[— (r — 80 AU) 2 /(2u; 2 )]} exp[— (r/k)]. (1) 

Since the primary heating source for the disk is stellar irradiation, the effects of shad- 
owing and illumination must be accounted for in determining the vertical structure of the 
disk. The three-dimensional density and temperature structure of the perturbed disk are 
calculated iteratively according to Jang-Condell (2008). That is, the illumination at the 
surface of the disk is determined by ray-tracing and the disk temperatures are calculated ac- 
cordingly. Once the disk temperatures arc determined, the vertical density profile of the disk 
is iteratively recalculated assuming hydrostatic equilibrium as above, keeping the vertically 
integrated surface density profile constant . 



4.3. Scattered Light 


The scattered light image of a disk is modeled as in Jang-Condell (2009) and Jang- 
Condell & Turner (2011). The scattering surface of the disk is defined to be where the 
optical depth from the star at a given frequency v is r v — 2/3. Then, the brightness at the 
scattering surface of the disk is 

rscatt _ 0 , v iiRlB v (%) 

4 r 2 (/i + cos i) 


where u) v is the albedo, p is the cosine of the angle of incidence to the scattering surface, B u 
is the stellar brightness, r is the total distance to the star, and i is the angle between the line 
of sight to the observer and the normal to the scattering surface. The observed brightness 
of a star at distance d is 
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so we can write express the surface brightness in scattered light in units of the apparent 
brightness of the star per square arcsecond: 
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Finally, one needs to choose a composition and size distribution of the dust to calculate 
the final scattered light that is emitted from the disk. Many compositions and possible grain 
sizes are available, and it is not immediately clear whether a single composition is preferred 
for optically thick disks, or whether each individual disk has its own properties. The neutral 
color of TW Hya in the optical compared to the quite red colors of HD 100546, suggest that 
a range of compositions and size distributions exist amongst disks. 

How predictive can models of optically thick disks be for composition and size distribu- 
tion? In an attempt to answer this question, we have constructed six independent structure 
models that vary composition between pure water ice and astronomical silicates (Warren 
1984; Laor & Draine 1993). and vary minimum grain sizes between 5x10 " 3 ,1, and 10pm. 
We take a fiducial disk model that has a maximum grain size of 1 cm, a gap width, w = 30 AU, 
gap depth, d— 0.3, located at a radius of 80 AU, and a truncation knee at fc=100 AU. To 
calculate the expected scattering efficiency of the dust grains as a function of wavelength, we 
calculated albedo and opacity using the program miex (Wolf & Voshchinnikov 2004). The 
dust is well-mixed with the gas so that the dust density is proportional to the gas density. 

Differing grain sizes and compositions can affect the resulting surface brightness profiles 
of a given disk structure. Figure 7 shows a comparison between the six predicted STIS surface 



brightness profiles for the different models. If the composition of the disk is incorrect, the 
structure of the disk can be incorrectly interpreted. 


In general, pure water ice disks are dimmer than those that possess pure silicates. 
The reason for this difference is primarily due to the higher opacity of a pure silicate disk, 
demonstrated in Figure 8. In this Figure we have plotted the height of the r„ = 2/3 surface 
( H ) divided by the radius ( R ) for varying wavelengths. From Eq. (4), the brightness of the 
disk is proportional to the angle of incidence n at the surface, 


JL ( H \ _ H _ r A. {—'] 

dR \ Rj ~ ~R ~ R Jr \r) 


(5) 


If H oc /? ', then fi = ( (3 — 1 )H/R. Therefore, the brightness scales roughly with H/R. The 
H surface occurs higher up in the disk for silicates, resulting in a. higher surface brightness. 
The different structure of the surface brightness profiles in Figure 7 is also striking, and is 
clarified by looking at d(H / R) / dR for the different models (See Figure 9). 


Water ice in particular has strong absorption features at 1.5 and 2.0 /wn both wave- 
lengths that are probed with our observations. Given the ~10% accuracy of the disk pho- 
tometry measurements, we can in principle detect absorption features in the medium band 
filters that have a depth of ~15-30%. Figure 10 shows the resulting reflectance spectra for 
each of our models. Silicates show no sharp features in the visible to Near-IR, but are mostly 
neutral over this wavelength range for minimum grain sizes >lp,m. Water ice shows notice- 
able features starting at ~l/im, but at the current levels of uncertainties in disk photometry, 
disks with minimum grain sizes of >10/nn would be required to detect such a feature. 


Based on our results above, composition can potentially play an important role in the 
observed scattered light spectrum as well as the observed surface brightness profiles of a 
disk. For our final model of TW Hya, an exhaustive search of compositions, grain sizes, and 
structures is computationally prohibitive and beyond the scope of this paper. We instead fo- 
cused on dust opacities and scattering efficiencies using a composition with the same relative 
abundances as D'Alessio et al. (2001): 29.69c organics, 40.4% water ice, 24.5% astronom- 
ical silicates, and 5.5% troilite. This model composition successfully matches the SED of 
TW Hya. We can then determine whether such a composition is warranted for the surface 
of the disk based on the scattered light observations. 


4.4. Model Parameters 


We ran a suite of models, varying the following disk parameters as follows: 

fl m in € (0.005. 0.5, 5} microns (6) 



w 


(7) 

( 8 ) 
(9) 
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G {10,20.30} AU 
d G {0.0, 0.3. 0.5} 
k G {60,80,100, 120,150} AU 

The surface brightness profiles are fit simultaneously spatially and across the 7 wave- 
lengths for which we have photometric data. To account for uncertainties in the overall 
normalization such as might arise from uncertainty in the distance to TW Hya, we allow the 
overall brightness to vary by a constant factor and calculate the reduced \ 2 value for each 
model with respect to the data from 56 to 160 AU. 

The best fit model is a m j n = 0.005 microns, w = 10 AU with a center at 80 AU, d = 0.3, 
and k = 150 AU. This has the lowest y 2 as compared to the observed data of 2.45. The 
model brightness profile must be scaled by an overall factor of 1.07 to fit the data, meaning 
that the data is 7% brighter than our models and can be explained by the uncertainty in 
TW Ilya’s distance which could account for brightness differences of up to 23%. 

Figures 11 and 12 show the radial data compared to the best fit model for all of the 
observations and radial spectral cuts of the disk. In general the model fits the observed 
spectrum quite well, with the exception of the region around 80 AU. This is also where the 
model fails to reproduce the observed surface brightness profile. 


5. Mass of the Gap-Opening Planet 


The best model for the TW Hya disk is one where an axisymmetric 30% partial gap is 
opened at 80 AU in the disk. A possible mechanism for opening such a gap is an embedded 
planet. The model we have used to fit to the TW Hya do not include hydrodynamics, but 
rather the gap is imposed in an ad hoc manner on the disk structure. 

Numerical hydrodynamic simulations of planets embedded in gas-dominated protoplan- 
etary disks by Bate et. al. (2003) indicate that tidal torques between the planet and disk can 
clear axisymmetric partial gaps. In the case of an inviscid disk, the gap-opening threshold 
is roughly the mass at which the Hill radius of the planet equals the thermal scale height of 
the disk. The Hill radius of the planet is 


nun 



( 10 ) 


where m p is the mass of the planet, M* is the mass of the star, and a is the star-planet 



( 11 ) 
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distance. The thermal scale height of the disk is 

rr C S 

H = — a 
y<p 

where c s = y/kT/fi is the local thermal sound speed, k is the Boltzmann constant, T is 
the local disk temperature, /x is the mean molecular weight, v<p — \J GM*/ a is the orbital 
speed of the planet, and G is the gravitational constant. For a composition dominated by 
molecular hydrogen, fi = 2rn H . 


For a viscous protoplanetary disk, a more massive planet is required to open a similarly 
sized gap in a disk. The relation can be expressed as 


3 H 50 

1 < l 

4 r Hiii qH ~ 


( 12 ) 


(Crida. et al. 2006), where q = m p /M * and 7 Z = r 2 kl P /u v is the Reynolds number, and u v is 
the viscosity, given by v v — ac s H for an o-disk model. 


Bate et al. (2003) adopt disk parameters of H/r — 0.05 and 1Z = 10 5 for their simu- 
lations, giving a gap-opening threshold of <7 = 1.06 x 10 ” 3 using Eq. (12), or slightly more 
than 1 M j. For comparison, the inviscid gap-opening threshold would be 0.4 Mj. They find 
that planets with q = 3 x 10 ~~ 4 and 1 x 10 -4 clear gaps of 90% and 50%, respectively, and a 
planet with q — 3 x 10 ~ 5 creates a nearly negligible gap. Thus, a planet that clears a gap of 
30% would be between 0.03 and 0.1 of the viscous gap-opening threshold. 


For TW Hya, our best-fit disk model has H/r = 0.077 and a = 2.5 x 10~ 4 , giving 
1Z = 7 x 10° and a viscous gap opening threshold of q = 7.8 x 10~ 4 . Since M* = 0.4 M e , 
this is 103 M®. This implies that if an embedded planet is the cause of the gap we have 
modeled, it must be between 3 — 10 M®. However, there is some degeneracy with regard 
to a, as discussed in §4.1. A larger value of a gives a lower Reynolds number, requiring a 
more massive planet to open a gap and allowing a larger planet to hide in the disk. If we 
adopt a = 5 x 10~ 4 , corresponding to a maximum grain size of 1 mm, we find a gap opening 
threshold of q = 10“ 3 . Thus, a more conservative upper limit on the mass of a planet at 80 
AU in the TW Hya disk is 13 M®. If the mass of TW Hya is larger, these numbers will need 
to be revised upwards. However, a larger stellar mass also makes our disk models flatter and 
dimmer, which does not fit our scattered light images. 


6 . Discussion 


We have combined several resolved images of the TW Hya disk in scattered light with 
spatially resolved spectroscopy of the disk. The spectrum of the disk in the visible and near- 
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1R is mostly featureless, with a broad neutral to blue trend, indicative of sub-micron grains. 
Based on our scattering models of the disk, a composition made primarily of organics, water 
ice, and silicates with grain sizes as small as ~0.005/«n best match the available scattered 
light data, and are consistent with models of TW Ilya’s SED. 

We can compare our predicted water ice mass abundance relative to the models of 
Hogerheijde et al. (2011) based on their detection with Herschel of water vapor at 100- 
150 AU in the TW Hya disk. They interpreted the deficit of water vapor to that predicted 
by their models as evidence for dust settling of large icy grains to the disk midplane and 
assumed that there needed to be > 9 x 10 27 g of water ice in the disk. Our scattered light 
images probe deeper in disk scale height than UV by definition and are more consistent with 
increasing water ice abundance, which would support the idea of dust grain settling. If we 
assume a disk mass in gas of 1.9xl0~ 2 M @ , assumed by Hogerheijde et al. (2011), then we 
would predict a total ice mass in the disk from our models of 2xl0 2 9 g, based on our mass 
abundance of 40% water ice in dust and a dust-to-gas ratio of 0.0138. Assuming that we 
probe the local mass abundance of water ice relative to gas at our r = 2/3 surface at the 
disk, this would imply a relative abundance of 5.6 xlO -3 at 30 AU above the midplane at 
100 AU and 46 AU above the midplane at 150 AU. 

Our models, which include a 30% cleared gap at 80 AU and a disk truncation exterior 
to 150 AU match the observed surface brightness profiles as a function of wavelength. We 
have placed an upper limit on the mass of a planetary companion that could be clearing 
the gap seen in TW Hya to between 5 — 13 M®. Despite a general agreement between our 
models and the data, we discuss here some of the discrepancies and their potential causes. 

The advantage of scattered light measurements is the ability to constrain minimum grain 
size since the smallest dust typically has the largest available surface area for scattering. Our 
models allow minimum grain sizes between 0.005 and 0.5 /mi with nearly equivalent chi-sq 
values. By the absence of any strong water ice absorption features, we can constrain the 
abundance of water ice only on large grains; small grains show absorption features at <5% 
the disk scattering continuum even if wholly made of ice (see Figure 10), below the level of 
our disk photometric uncertainties. We have chosen a composition based on what is most 
likely present in the TW Hya disk, but we have not exhaustively explored other materials. 
Any material that scatters neutrally in the visible to near-IR would most likely also be 
suitable to explain the spectral shape we observe. To test multiple grain sizes requires a 
self-consistent model with different parameters used by D'Alessio et al. (2001). Our current 
models for TW Hya are being expanded to include these effects to better model the disk 
(Jang-Condell 2011, in prep). 

Besides color, the other main effect of grain size is on the overall surface brightness of the 



disk, and the higher albedo of smaller grains helps match the data. Yet, our models slightly 
underpredict the measured brightness; we need to scale the model by a factor of ~1.07. The 
presence of higher opacity grains of neutral color, possibly the presence of more silicates 
in the small grains, could make the disk brighter yet. This may argue for fewer organics 
and ices in the TW Hya disk. Dust of size comparable to the wavelength is expected to 
forward-scatter; however, including forward scattering would decrease the overall brightness 
of the disk because we observe mainly 90° scattering. 

Interestingly, the brightness asymmetry between the front, and back sides of the disks is 
only apparent at wavelengths <l.l/nn, and the brightness asymmetry for wavelengths longer 
than 1.1/iin are not well fit by our models if we assume some disk inclination. The short 
wavelengths probe the uppermost disk surface where the bulk of the vertical height of the 
disk is optically thick to them. This suggests that the brightness asymmetry is caused by a 
perturbation in the surface grains only. 

While our models are of the outer disk, which is mainly illuminated by the central star 
and not affected by accretion, the general framework of the a viscosity disk does set the 
disk viscosity and therefore the ease with which a gap is opened. The upper limit of the 
gap-opening mass depends approximately linearly on the viscosity, which in turn is derived 
from the measured accretion rate. 

Finally, we investigate other possibilities for the presence of the gap. It is interesting 
to note in Figure 12 that at 80 AU, the center of the gap, the spectral character of the 
disk is different than that seen in other portions of the disk. If this is real, then a different 
composition in this annulus may be causing the overall structure of the disk to change. 
Conversely, compositional changes in this gap might be caused by enhanced collisions within 
the disk or the presence of an accreting protoplanet. Both of these possibilities are supported 
in part by the recent discovery to a sharp cutoff in the millimeter emission from dust at 
60 AU, which is not present in similar high resolution imaging of CO gas (Andrews et al. 
2011). Since we see scattered light from dust well beyond this cut-off it is clearly due to a 
change in the bulk size distribution or millimeter emission properties of the dust, which in 
of itself may be a signature of a planetary companion or change in bulk dust composition. 

The gap could instead be an unresolved spiral structure in the disk, as seen in HD 141569 
(Clampin et al. 2003). This could explain both the asymmetry as well as the decrease in 
surface brightness. Higher spatial resolution images, such as from next generation ground 
based AO coronagraphs, might solve the mystery. But. the rather small field of view for these 
instruments (~2"= 112 AU) might not be able to probe these interesting features in the TW 
Hya disk effectively. 
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F204M F222M 


Fig. 1. - Final, roll combined. PSF-subtractccl images of TW Hya in the four medium band 
filters observed, FI 71 M. F180M, F204M, and F222M. Regions in which there is no data, 
i.e. under the coronagraph plus in an extended region in the NE/'SW diffraction spike are 
blacked out. The mottled appearance in the shorter filters, particularly F204M, is the result 
of systematic residuals from diffraction and scattering patterns. The bright diagonal stripe 
at lower SW is from the edge of the NIC-2 array. 
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Fig. 2. — Final images from previous scattered light observations of TW Hya with ST1S, 
NICMOS’ FI 10W filter, and F160W filter. 
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Fig, 4. Azimuthal brightness profiles for the STIS and F1KJW images between 50-130 AU 
(top) and 130-230 AU (bottom). Each curve was normalized by the median surface bright- 
ness. In each panel we compare the brightness asymmetry with that of our best fitting model 
disk in §4.4 with inclinations of 15° (solid line) and 7° (dashed line). 
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Fig. 5. Azimuthal brightness profiles for the F160W and F171M images between 50- 
130 AU (top) and the F180M and F222M images between 50-130 AU (bottom). Each curve 
was normalized by the median surface brightness. In each panel we compare the brightness 
asymmetry with that of our best fitting model disk in §4.4 with inclinations of 15° (solid 
line) and 7° (dashed line). 





star 



0.007 


0.006 


0.005 


0.004 


0.003 



0.5 1.0 1.5 2.0 2.5 

Wavelength (jum) 


Fig. 6. - Total spectral reflectance of the TW Hya disk from 0.5-2.22/nn, including an 
average of the spatially resolved spectra of TW Hya from Roberge et al. (2005) normalized 
to the STIS disk photometry. 
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Fig. 7. ----- Comparison of TW Hya-like disk surface brightness profiles for different pure water 
ice (squares) and pure astronomical silicate (solid and dashed lines) compositions at, the STIS 
bandpass. 
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Fig. 8. A comparison of H/'R for astronomical silicates (blue lines) and water ice (red lines) 
with a m j n ==0.005//.m. where H is the r = 2/3 optical depth surface at each wavelength. As 
expected, longer wavelength filters probe deeper within the disk surface, though the specific 
geometry of that surface changes with composition. 
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Fig. 9. - Comparison of d(H/R)/dR for pure water ice and astronomical silicate composi- 
tions of differing a min . The derivative of the disk surface goes into the calculation of scattered 


light as discussed in §4.3. 
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Fig. 10. Comparison of TW Hya-like disk reflectance spectra for different pure water ice 
(squares and asterisks) and pure astronomical silicates (solid and dashed lines). Squares 
and solid lines represent the interior of the disk at 56 AU, while asterisks and dashed lines 
represent the outer part of the disk at 140 AU, 
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Fig. 12. Comparison of radial spectra of the disk with model profiles from the best fitting 

disk model of §4.4. The spectra are normalized to the brightness at the STIS wavelength of 
0.58 pm and offset vertically by whole numbers, starting with 3 for 60 AU, 2 for 80 AU, 1 
for 120 AU, and 0 for 140 AU. 
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Table 1 . Photometry of TW Hya and PSF Scalings 


Filter 

PSF Reference Star 

Scaling 

TW Ilya Flux Density (Jy) 

STIS 50CCI) 

HD 85512 

0.0544,0.0586 

0.1955 

FI 10W 

r 1 Eridani 

0.099T0.0027 

0.747 

F160W 

G1 879 

0.031T0.002 

1.03 

F171M 

CD-43°2742 

0.541T0.008 

0.94T0.01 

F180M 

CD-43°2742 

0.528±0.001 

0.84T0.01 

F204M 

CD-43°2742 

0.56T0.01 

0.77T0.01 

F222M 

CD-43°2742 

0.61T0.01 

0.76T0.02 



